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Abstract. The 6-12 /xm spectrum of Neptune has been 
recorded with the PHT-S instrument of the Infrared Space 
Observatory (ISO) at a resolution of 0.095 /xm. In addi- 
tion to the emissions of CH 4 , CH 3 D and C 2 H 6 previously 
identified, the spectrum shows the first firm identifica- 
tion of ethylene C 2 H 4 . The inferred column density above 
the 0.2-mbar level is in the range (1.1-3) xlO 14 molecules 
cm” 2 . To produce this low amount, previous photochem- 
ical models invoked rapid mixing between the source and 
sink regions of C 2 H 4 . We show that this requirement can 
be relaxed if recent laboratory measurements of CH 4 pho- 
tolysis branching ratios at Lyman a are used. 
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1. Introduction 

In the stratospheres of all giant planets, methane (CH 4 ) 
photochemistry leads to the formation of many hydrocar- 
bons. Photochemical models have been developed by sev- 
eral authors (e.g. Gladstone et al. 1996 for Jupiter, Moses 
et al. 1999 for Saturn, Atreya et al. 1991 for Uranus, and 
Romani et al. 1993 for Neptune) to account for the obser- 
vations of the detected hydrocarbons (C 2 H 2 and C 2 H 6 ) 
and to predict the abundances of minor species. Among 
the expected products in significant amount is ethylene 
(C 2 H 4 ). This compound however was only detected in the 
northern auroral region of Jupiter, from emission in the ^7 
band at 10.5 /xm observed in Voyager/IRIS spectra (Kim 
et al. 1985). On Neptune, a tentative detection of the same 
band of ethylene was presented by Orton et al. (1987) from 
ground-based low-resolution spectral measurements. 

Observations with the short- wavelength spectrometer 
(SWS) of the Infrared Space Observatory (ISO) have led 
to the detection of several new hydrocarbons in Saturn 
(C 3 H 4 , C 4 H 2 , and CH 3 ; de Graauw et al. 1997; Bezard et 
al. 1998) and Neptune (CH 3 ; Bezard et al. 1999b). In the 
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case of Neptune, the 7-14 /xm SWS spectrum, recorded at a 
spectral resolution of 0.005 /xm has been studied by Bezard 
(1998). Information was retrieved upon the stratospheric 
temperature, the CH 4 stratospheric mixing ratio, and the 
C 2 H 2 and C 2 H 6 abundance profiles. In the case of C 2 H 4 , 
only an upper limit of its stratospheric column abundance 
( 8 xl 0 14 molecules cm ” 2 above the 0 . 2 -mbar level) was 
inferred (Bezard et al. 1999a). 

We present in this paper the 6-12 /xm spectrum of Nep- 
tune recorded with another ISO instrument, ISOPHOT, 
and we report the first firm detection of C 2 H 4 on Nep- 
tune. Section 2 describes the observations and the data 
reduction. Section 3 presents the modelling and the deter- 
mination of the C 2 H 4 abundance. Results are discussed in 
Section 4 and conclusions are presented in Section 5. 

2. Observations and data reduction 

The photometer ISOPHOT (Lemke et al. 1996) of ISO 
(Kessler et al. 1996) included a subsystem (PHT-S) made 
of two grating spectrometers operating at 2.5-5/xm (SS) 
and 6-12 /xm (SL) with a spectral resolution of 0.0445 /xm 
and 0.0949 /xm respectively (Klaas et al. 1997). The two 
spectra were simultaneously recorded on the two linear 
arrays of 64 elements Si:Ga detectors each. The entrance 
aperture was 24 ;/ x24". 

The PHT-S spectrum of Neptune was recorded on May 
8 , 1997, between 3:41 UT and 12:29 UT. The diameter of 
Neptune was then 2.25", and the sub-Earth latitude was 
-30.6°. Twelve pairs of spectra, of 1024 sec integration 
time each, were successively recorded on and off the source 
(with the offset position at 5' north), in order to subtract 
the sky background. 

The data reduction was performed using the PHT 
Interactive Analysis software package PI A V7.3.2(e) 
(Gabriel et al. 1997) together with IDL routines that facil- 
itate individual inspection of the signals and deglitching. 
The Si:Ga photoconductors, arranged in two arrays of 64 
pixels each, exhibit transients of the photocurrent after a 
flux change. This behavior impacts the deglitching proce- 
dure as well as the assignment of flux densities. 
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Fig. 1 . PHT-S spectrum of Neptune observed in the 6-11.6 /jm 
range. Emission features from CH4, CH3D, C2H6 and C2H4 
are dearly detected. The flux scale on the left axis corresponds 
to the direct PHT-S calibration, that on the right axis to a 
rescaling of the data by a factor of 0.67 to agree with the SWS 
observations of the 7.7-/zm CH4 band. 


The 32 s integration ramps with 128 readouts were 
divided into 16 subramps. The fitted slopes s(f), being 
a measure of the photo current, were deglitched auto- 
matically using an algorithm that clips all values where 
| $(t) — £>(<)| > 3.5<r. b(t) are the boxcar smoothed slopes 
with the smoothing parameter tuned to follow the slow 
transients of the photo current w ith time t. a is calcu- 
lated as — b(t)) 2 /(n — 1), where n is the number 

of slopes. 

The assignment of flux densities was performed us- 
ing the new dynamic calibration method (Acosta-Pulido 
1999). It is based on the reproducibility of the transient 
behavior when a given detector pixel is exposed to sim- 
ilar flux levels. The method, which is implemented into 
PI A, requires the prior subtraction of the exposure time- 
dependent dark current. The transients also differ from 
pixel to pixel, so PIA contains a database for all pixels of 
measured transients with 27 calibration sources, covering 
most of the accessible flux range. After deriving 12 spectra 
on source and 12 on the background, we concatenated the 
pairwise differences within PIA and then averaged each 
pixel to derive the final flux density. 

The uncertainties in the PIA results and propagated 
statistical errors are usually smaller than the variations 
due to long-term responsivity changes and transient re- 
peatability. Consistency checks within the database of 
transients suggest average accuracies of 10 to 15% over all 
pixels. This is also consistent with results from long-term 
stability checks (Schulz 1999) which give 10%. Relative un- 
certainties between fluxes of adjacent detector pixels with 
similar flux levels are smaller, i.e, closer to the statisti- 
cal errors. The error in the SED models of the calibration 
standards is not included and ranges between 3 and 5%. 


Fig. 1 shows the PHT-S spectrum of Neptune between 
6 and 12 /zm. The 1-cr noise level is estimated to about 
0.010 Jy at 6-7 /im and 0.020 Jy at 10-11 /im. The spec- 
trum is dominated by the strong z/ 4 band of CH 4 which 
peaks at 7.7 /im. CH4 emission is also seen in the Q-branch 
of the weak v 2 band at 6.5 /im, observed here for the first 
time. CH3D is detected at 8.6 /im through its 1 /$ band 
and contributes with its band to the emission at 6.8 
/im. The wing of the z/9 band of C2H6 is visible beyond 

11.4 /im. All these emissions above 7 /im were previously 
detected in the SWS spectrum of Neptune (Bezard 1998). 
In addition, the PHT-S spectrum of Neptune shows a new 
emission feature in a single pixel, at 10.5 /zm, at a level of 
0.1 Jy, i.e. 5 times above the noise level. We attribute this 

10.5 /zm emission to the Q-branch of the C2H4 v 7 band, 
which is centered at 949 cm -1 (10.53 /zm). 

The comparison of the PHT-S and SWS data shows a 
very good agreement in the shapes of the two spectra, but 
a difference in the absolute calibration scale by a factor 
of 1.5, the PHT-S fluxes being higher than the SWS ones. 
The two scales are indicated on Fig. 1. We have presently 
no explanation for this discrepancy. 

In this paper, we modelled the observations with the 
“nominal” atmospheric model (methane and temperature 
profiles) derived by Bezard (1998) from an analysis of the 
SWS spectrum between 7 and 14 /zm (see Section 3). The 
C2H2 and C2H6 abundances inferred from this model were 
in very good agreement with previous results from Voyager 
data (Bezard et al. 1991) and ground-based infrared mea- 
surements (Orton et al. 1987). For consistency, we used 
here the SWS calibration of the PHT-S data to infer the 
“nominal” C2H4 abundance, but we took into account the 
difference in the two instrument calibrations in the error 
analysis. 

3. Radiative transfer analysis 

Synthetic spectra were generated from a line-by-line radia- 
tive transfer program that includes the H2-He collision- 
induced absorption and molecular opacity from CH4, 
CH3D, C2H6, and C2H4. Line parameters come from the 
1997 spectroscopic Geisa databank (Jacquinet-Husson et 
al. 1997). We used the same temperature profile as in 
Bezard (1998) and Bezard et al. (1999b). C2H4 spectral 
calculations were carried out assuming LTE at all atmo- 
spheric levels. 

To derive constraints on the C2H4 stratospheric col- 
umn abundance, we proceeded as did Bezard et al. in their 
analysis of the methyl radical emission. We considered ver- 
tical profiles derived from photochemical calculations, al- 
lowing for the uncertainty in the photochemical scheme 
(Profiles 1-5 in Table 1) and in the eddy mixing pro- 
file (Profiles 6-7) (see Section 4). To determine the C2H4 
abundance implied by the observations, we then multi- 
plied each test profile by a constant factor that allowed us 
to best reproduce the intensity of the observed 10.5-/zm 
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emission in a least-squares sense. We found column den- 
sities for these rescaled profiles in the range 1.7-2.2xl0 14 
molecules cm -2 down to the 0.2-mbar level (Table 1). The 
C2H4 abundance below this level does not contribute sig- 
nificantly to the emission, being located at colder levels in 
the stratosphere. 

The spectrum generated with our nominal model (Pro- 
file 5 for C2H4) is shown in Fig. 1 in the whole range 6-11.6 
fim. The discrepancy between the synthetic and observed 
spectra around 6.8 fim may be due to the presence of the 
v$ band of C2H6. This band was not included in the calcu- 
lations as no line compilation exists in the usual spectral 
databanks. In the other regions, the agreement is quite sat- 
isfactory. Spectra incorporating various vertical profiles of 
C2H4 (Profiles 1-5, Fig. 2) are shown in Fig. 3 (see Section 
4 for a full description of the profiles). 

Our nominal value for the inferred column density is 
1.7 x 10 14 molecules cm -2 , using the rescaled profile 5. Un- 
certainties arise from the noise level (±20%), the flux cal- 
ibration of the ISO spectrum, the model temperature pro- 
file, and the C2H4 vertical profile. For the flux calibration 
uncertainty, we assumed a ^30% error interval that in- 
cludes the PHT-S flux scale. Uncertainties from the tem- 
perature profile were estimated using the extreme “cool” 
and “warm” profiles from Bezard et al. (1999b). A l 2 g% 
error range was then obtained. Finally the uncertainty 
due to the C2H4 height profile is ig 4 5 * % (Table 1). Com- 
bining the above error bars, we conclude that the C2H4 
column density above the 0.2-mbar level is in the range 
(1.1-3) xlO 14 molecules cm" 2 . 

4. Comparison with photochemical models 

Photochemical calculations were carried out with a one- 
dimensional model as in Bezard et al. (1999b). The model 
has been most recently described in Bishop et al. (1997). 
Chemical reactions and kinetic rates included in the mod- 
eling are listed in Table 3 of Bishop et al. with minor 
updates since then, except for the methane photolysis 
branching ratios at Lyman a and the rate coefficient for 
C2H4 ± H ± M C2H5 ± M, which are both discussed 
below. Photolysis rates are calculated for disk-averaged 
conditions, and account for both solar irradiance and the 
Lyman a skyglow from the local interstellar medium. So- 

lar minimum conditions, representative of the time of the 
ISO observations, were used, because the chemical lifetime 
of C2H4 was found to be shorter than the solar cycle. We 
used a CH4 mole fraction equal to 1.4xl0 -3 in the lower 
stratosphere as favored by observations of the v\ band of 
methane by ISO (Bezard 1998). The model incorporates a 
downward flux of atomic hydrogen at the upper boundary 
equal to 4x 10 7 cm" 2 s _1 , representing H production from 
solar EUV at solar minimum. 

We used as “nominal” an eddy mixing coefficient, K, 
proportional to the atmospheric number density n to the 
-0.6 power, and equal to K/j = 10 7 cm 2 s -1 at the methane 
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Fig. 2 . Temperature-pressure profile used in the radiative 
transfer analysis (solid line with dots, top axis) and calculated 
C 2 H 4 profiles. Three different CH 4 photolysis branching ratios 
at Lyman a were first used: Romani (1996) pong-dashed line], 
Heck et al. (1996) [dash-dot line], Smith and Raulin (1999) 
[solid line], all with the nominal value of the reaction rate co- 
efficient for H + C 2 H 4 from Lightfoot and Pilling (1987). Dif- 
ferent expressions of the reaction rate coefficient of H ± C 2 H 4 
were also considered: Lightfoot and Pilling’s (1987) nominal 
expression [solid line], Lightfoot and Pilling’s fast rate [short- 
dashed line], Baulch et al. (1994) [dotted line], all with CH 4 
photolysis branching ratios from Smith and Raulin (1999). 

homopause. This type of K profile was shown to be com- 
patible with the ISO observations of CH3 if the low pres- 
sure rate coefficient for the methyl self-reaction was be- 
tween that of Slagle et al. (1988) and Macpherson et al. 
(1985). Following Bezard et al. (1999b), we also varied the 
K profile, considering constraints from other observations 
of Neptune hydrocarbons (Bishop et al. 1997). 

Ethylene photochemistry is similar to that of the 
methyl radical CH3 (see Bezard et al. 1999b). It has a 
source high in the stratosphere (3 to 0.3 /ibar region) 
from methane photolysis, a three-body sink low in the 
stratosphere (0.1 to 1 mbar region), and transport con- 
necting the two. Methylidyne (CH) is directly produced 
from methane photolysis as is methyl, and yields ethylene 
via CH ± CH4. Triplet methylene ( 3 CH2), another prod- 
uct of methane photolysis, also leads to C2H4 formation. 
The net effect is that C2H4 production is sensitive to the 
sum of the photolysis pathways that yield either 3 CH2 or 
CH. The dominant sink for ethylene is, in the lower strato- 
sphere, the three-body reaction C2H4 ± H + M — ► C2H5 
± M . The C2H4 column density is proportional to its 
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* production rate but inversely proportional to the loss rate 
coefficient and to the transport rate between the source 
and the sink. 

The methane photolysis rate is dominated by the flux 
at Lyman a. Unfortunately, the branching ratios for the 
possible pathways at Lyman a are still not known. Fur- 
thermore, there are disagreements among the laboratory 
experiments about the relative importance of the different 
pathways and the quantum yield of H and H 2 in methane 
photolysis (Smith and Raulin 1999, Brownsword et al. 
1997, Heck et al. 1996, Mordaunt et al. 1993). We first 
used the same branching ratios as in Romani (1996) (Pro- 
file 1 in Fig. 2). These branching ratios are based upon 
the work of Mordaunt et al. (1993) and previous work. As 
can be seen in Fig. 3, this produces too large a column 
density of C 2 H 4 to be compatible with the observations. 
We next tried the branching ratios suggested by Heck et 
al. (1996) (Profile 2). The main difference with those of 
Romani (1996) is that the production of 3 CH 2 + 2H is 
reduced from 0.21 to zero while the production of CH 3 + 
H is increased by a similar amount. The column density of 
C 2 H 4 decreases, but not enough to match the observations 
(see Table 1 and Fig. 3). Lastly we used the branching ra- 
tios of Smith and Raulin (1999) (Profile 3), which are very 
similar to those of Brownsword et al. (1997). With respect 
to ethylene formation, the major difference between these 
branching ratios and those of Heck et al. (1996) is a re- 
duction in the branching ratio to form CH + H 2 + H from 
0.11 to 0.06. The modelled C 2 H 4 column density is now 
close to the value derived from the observations. 

The predicted C 2 H 4 emission is still slightly too high, 
although the branching ratios of Smith and Raulin re- 
sult in the lowest model predicted C 2 H 4 column density 
of all that we tried. We next investigated the loss pro- 
cesses for ethylene. We used the “fast” rate coefficient for 
H H- C 2 H 4 from Lightfoot and Pilling (1987) (the largest 
possible rate coefficient within their reported error bars) 
(Profile 4). Using this rate coefficient with the Smith and 
Raulin branching ratios, the modelled C 2 H 4 profile now 
yields the correct abundance to reproduce the observed 
emission within the noise level (see Table 1 , Figs. 2-3). 
It can be noted that the rate coefficient expression sug- 
gested by Baulch et al. (1994) for this reaction is faster 
than the “fast” Lightfoot and Pilling (1987) value. It re- 
sults in a predicted C 2 H 4 emission still consistent with the 
observations (Profile 5). 

It is also possible to reduce the C 2 H 4 column density 
by changing K. A first way is to keep the same func- 
tional form, K~n _0 6 , but to decrease its value at the 
methane homopause. Methane photolysis then occurs at 
higher number densities, closer to the C 2 H 4 loss region. 
To reproduce the C 2 H 4 emission feature with the Smith 
and Raulin (1999) branching ratios K must be lowered to 
5xl0 6 cmV 1 at the methane homopause, but with the 
Romani (1996) branching ratios K must be decreased be- 
low 10 6 cmV 1 . The former value is within the bounds 



Fig. 3. Synthetic spectra calculated with the C 2 H 4 profiles 1-5 
of Table 1 (solid lines), displayed in Fig. 2 , are compared with 
the PHT-S spectrum in the region of the 1^7 ethylene band 
(histograms). 

permitted by the ISO/S WS CH 3 observations while the 
latter is not (Bezard et al. 1999b). A second way is to 
alter the form of K to increase the downward transport 
of C 2 H 4 from its source to sink. This was done previ- 
ously by Romani et al. (1993) and Bishop et al. (1997). 
They found that the model predicted C 2 H 4 with K~n " 0,6 
was too high for the Voyager UVIS solar occultation light 
curves (note that these papers were prior to Brownsword 
et al. (1997) and thus had high C 2 H 4 production). They 
solved the problem by increasing the transport rate be- 
tween the /xbar and mbar region (type “B” eddy profile in 
Romani et al. 1993) and using the fast rate coefficient for 
H -|- C 2 H 4 from Lightfoot and Pilling. With this combina- 
tion of transport and chemistry, the C 2 H 4 column density 
is now within the range allowed by ISOPHOT although 
at the upper bounds. 

5. Summary and conclusions 

The PHT-S detection of the C 2 H 4 j /7 band on Neptune 
implies a column density equal to 1.71 q 3 x10 14 molecules 
cm” 2 . This result is fully consistent with the upper limit of 
8xl0 14 molecules cm -2 reported by Bezard et al. (1999a) 
from ISO-SWS observations. Orton et aids (1987) obser- 
vations of Neptune at a resolution of 0.23 /im showed en- 
hanced emission around 10.5 //m that was attributed to 
C 2 H 4 . Their reported brightness temperature at 10.5 fim 
corresponds to a 0.06 Jy flux (±20%), in good agreement 
with the PHT-S spectrum which yields a value of 0.09 
(PHT calibration) or 0.06 Jy (SWS calibration) after con- 
volution at 0.23-//m resolution. 

For the photochemical model to reproduce both the 
ISO CH 3 and C 2 H 4 emission features with K~n “ 0 6 , the 
CH 4 photolysis branching ratios at Lyman a must be 
similar to those proposed by Brownsword et al. (1997) 
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Table 1 . Comparison of photochemical models to ISO observations 


Model K h 

(cm 2 s -1 ) 

CH 4 branching ratios 
at Ly a 

H -f C 2 H 4 rate coefficient 

Model C 2 H 4 Scaled C 2 H 4 a 
(molecules cm -2 at 0.2 mbar) 

1 

lxlO 7 

Romani (1996) 

Lightfoot k Pilling (1987) nominal 

8.8xl0 14 

1.9xl0 14 

2 

lxlO 7 

Heck et al. (1996) 

Lightfoot k Pilling (1987) nominal 

4.3xl0 14 

1.9xl0 14 

3 

lxlO 7 

Smith k Raulin (1999) 

Lightfoot k Pilling (1987) nominal 

3.0x 10 14 

1.9x 10 14 

4 

lxlO 7 

Smith k Raulin (1999) 

Lightfoot k Pilling (1987) fast 

2.2xl0 14 

1.8xl0 14 

5 

lxlO 7 

Smith k Raulin (1999) 

Baulch et al. (1994) 

1.6xl0 14 

1.7 x 10 14 

6 

1x10® 

Smith k Raulin (1999) 

Lightfoot k Pilling (1987) nominal 

2.2x 10 14 

2.2xl0 14 

7 

5xl0 7 

Smith k Raulin (1999) 

Lightfoot k Pilling (1987) nominal 

3.3xl0 14 

1.8xl0 14 


a Needed to reproduce the C 2 H 4 emission intensity 


and Smith and Raulin (1999). Furthermore, if K* = 10 7 
cmV 1 , the rate coefficient for H + C 2 H 4 must be near 
the “fast” value of Lightfoot and Pilling (1987). 

It is also possible to invoke rapid transport between the 
//bar and mbar region (type “B” K profile from Romani 
et al. 1993). In this case, the CH 4 photolysis branching 
ratios at Lyman a that lead to C 2 H 4 production must 
be increased above those proposed by Brownsword et al. 
(1997) and Smith and Raulin (1999). However, a fast rate 
for H + C 2 H 4 is still needed. 

In summary, previous photochemical models required 
rapid transport and a fast sink for C 2 H 4 to reproduce ob- 
servations. Now the new laboratory measurements of CH 4 
photolysis branching ratios at Lyman a allow us to relax 
that requirement. The need for laboratory measurements 
mentioned by Romani et al. (1993) — the CH 4 branching 
ratios at Lyman a, and the rate coefficient of H + C 2 H 4 
at low temperatures and pressures — is still present. 
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